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Abstract We revisit the well-known Hoyle-Narlikar- 
Whceler proposition that neutron star emerging in the 
magnetie-flux-conserving proeess of eore-collapse su- 
pernova can eonvert the stored energy of Alfven vibra- 
tions into power of magneto-dipole radiation. We show 
that the necessary requirement for the energy conver- 
sion is the decay of internal magnetic field. In this case 
the loss of vibration energy of the star causes its vibra- 
tion period, equal to period of pulsating emission, to 
lengthen at a rate proportional to the rate of magnetic 
field decay. These prediction of the model of vibra- 
tion powered neutron star are discussed in juxtaposition 
with data on pulsating emission of magnetars whose ra- 
diative activity is generally associated with the decay 
of ultra strong magnetic field. 

Keywords neutron stars, torsion Alfven vibrations, 
magneto-dipole radiation 



1 Introduction 

It is generally known today that the existence of neu- 
tron stars with dipole magnetic field of up to i? = 
10^'* — 10^^ G, emerging in the magnetie-flux-conserving 
core-collapse supernova, has first been regarded by 
Woltjer (1964). Based on this Hoyle, Narlikar and 

Sergey Bastrukov 
Junwei Yu 
Irina Molodtsova 
Renxin Xu 

State Key Laboratory of Nuclear Physics and Technology, School 
of Physics, Peking University, Beijing 100871, China 

Joint Institute for Nuclear Research, 141980 Dubna, Russia 

^ State Key Laboratory of Nuclear Physics and Technology, 
School of Physics, Peking University, Beijing 100871, China 

^ Joint Institute for Nuclear Research. 141980 Dubna, Russia 



Wheeler (1964) proposed that a neutron star with such 
a field can generate magneto-dipole radiation by means 
of conversion of energy of hydromagnctic, Alfven, vi- 
brations into power of electromagnetic emission (e.g., 
Pacini 2008). In recent work (Bastrukov et al. 2011) 
such a possibility has been studied in some details in the 
context of electromagnetic activity of magnetic white 
dwarfs. In this article, we present an extensive anal- 
ysis of vibration-energy powered magneto-dipole radi- 
ation of a neutron star undergoing global Alfven tor- 
sion seismic vibrations about axis of its dipole mag- 
netic moment. In §2, the general theory of node-free 
torsional Alfven vibrations of perfectly conducting de- 
generate stellar matter with frozen-in eonstant-in-time 
magnetic field is briefly outlined. In §3, the theory is 
extended to the ease of Lorentz-force-driven torsion vi- 
brations in time-varying magnetic field. It is shown 
that magnetic field decay is crucial (necessary condi- 
tion) to the energy conversion from Alfven vibrations 
to magneto-dipole radiation. The predictions of the 
model of vibration powered neutron star are discussed 
in juxtaposition with available data . The results are 
summarized in §4 with emphasis on the relevance of the 
model to electromagnetic activity of magnetars. 



2 Alfven node-free vibrations in 

constant-in-time frozen-in magnetic field 

The starting point in the study of node-free Alfven 
(non-compression magneto-mechanical) vibrations of a 
neutron star (thought of as a spherical solid mass 
of a perfectly conducting non-fiowing continuous mat- 
ter with frozen-in eonstant-in-time magnetic field) are 
equations of magneto-solid-mechanics (Bastrukov et al. 
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2009a, 2009b, Molodtsova et al. 2010) 
pu = -[Sj X B], 



(1) 



(5j = — [V X (5B], (5B = Vx[uxB], (2) 
An 

pii = -^-[V X [V x [u x B]]] X B, V • u = 0. (3) 

47r 

In such vibrations the rate of differentiaUy rotational 
material displacements of stellar matter is described by 
nodcless toroidal vector field 

u(r, t) = [a;(r, t) x r], uj{r, t) = [Vx(r)] a{t), (4) 

(5) 
(6) 



V\{r) = 0, xir)=feir)Piicos9), 



which is identical to that for torsion node-free vibra- 
tions restored by Hooke's elastic force (Bastrukov et al. 
2007a, 2007b). In the last equation, P^(cos6') stands 
for Legendre polynomial of degree i specifying the over- 
tone of toroidal a- mode; the amplitude a{t) describes 
time evolution of vibrations (both global and locked 
in the crust). In computing discrete spectrum of such 
vibrations, the magnetic field can be conveniently rep- 
resented in the form: B(r) = _Bb(r), where B is the 
intensity and b(r) is dimensionless vector- function of 
magnetic field distribution over the star volume. Sim- 
ilar representation can be used for the bulk density 
p{r) = p4>{r), where p is the density at the star center 
and 4>{r) describes the radial profile of density which 
can be taken from computations of neutron star struc- 
ture relying on realistic equations of state accounting 
for non-uniform mass distribution in the star interior 
(e.g., Weber 1999). 

Scalar product of ([3]) with the following separable 
representation of u{r,t) — a(r) a(t) and integration 
over the star volume leads to equation for a(i) having 
the form of equation of harmonic oscillator 

Ma(t)+K„,OL(t) = 0, 7W = pmg, /C^ = -;-fcf,(7) 



= y (/)(r) a(r) • a(r) dV, 



ki = 



J a(r) • [b(r) x [V x [V x [a(r) x b(r)]]]] dV, 



M 

ki 

it, VA 

mi 



R 



R 

3M' 



B 



The amplitude of oscillations with constant in time fre- 
quency is given by 



a{t) 



ao cos ojt, 



2Ea 2Ea 



Ea = {l/2)Ma^ 



(8) 
(9) 



where over-bar stands for averaging over period of vi- 
brations: a'^{t) = (l/2)aQ. This suggests, if all the en- 
ergy E'burst of X-ray outburst goes in the quake-induced 
vibrations and, i.e. when E'burst = Ea, the amplitude 
ao can be extracted from the last equation. In the re- 
minder of the paper we remove index i and confine our 
analysis to the case of quadrupole overtone of a-mode, 
i.e., putting uj = Wf=2- The important outcome of as- 
sumption about constant in time undisturbed magnetic 
field is the vibration energy conservation 



dEAjt) 
dt 



0, EA{t) 



(10) 



and also that the fundamental frequency j/^ — ^jaI^-k 
(where w^i — va/ R) and the period Pa 
Alfven oscillations 



Vji^ ^ of global 



B / R 



27r /3M 



(11) 



remain constant in time. The magnitudes of basic fre- 
quency of toroidal a mode in neutron stars with mag- 
netic fields typical to radio-pulsars B12 = i3/(10^^ G 
and magnetars (soft gamma repeaters) i?i4 = B/{1Q^^ 
G are 

va{Bi2) = 2.055 • 10-3Bi2i?6/'(Af/Afo)-i/2^ (12) 
va{Bia) = O.2O55-Bi4i?6^^(A-^/M0)"^/2 (^3) 




Fig. 1 The magnetic field decay and period elongation 
in a quaking neutron star undergoing Lorentz-force-driven 
quadrupole torsion vibrations caused by conversion of vi- 
bration energy into energy of magneto-dipole emission. 
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where Re = i?/(10^ cm. It should be noted that above 
equations of sohd-star model are applicable not only to 
neutron stars but also for white dwarfs (Molodtsova et 
al. 2010; Bastrukov et al. 2010) and quark stars (Xu 
2003, 2009). In this work, continuing investigations re- 
ported in (Bastrukov et al. 2009a, 2009b, 2010), we re- 
lax the assumption about constant-in-timc undisturbed 
magnetic field and examine the impact of its decay on 
quake-induced vibrations and resultant magneto-dipole 
radiation. 



3 Energy conversion from Alfven vibrations 
into magneto-dipole radiation 

In the model under consideration, the decay of mag- 
netic field (during the time of post-quake vibrational 
relaxation of the star) is though of as caused by cou- 
pling of the vibrating star with outgoing material ex- 
pelled by quake. With allow for dependence on time 
of undisturbed magnetic field intensity, B = B(t), 
the above chain of argument leads to equation of vi- 
brations with the spring constant depending on time: 
^Aa{t) + K.m{B{t))a{t) = 0. The main subject of our 
further analysis is the impact of magnetic field decay on 
radiative process of neutron star undergoing Lorentz- 
forcc-driven vibrations obeying the last equation. It 
is easy to see that in this case the energy of vibra- 
tions is not conserved. The loss of vibration energy 
Ea = {l/2)[Md'^{t) + IC,n{B{t))a^{t) is determined by 
the rate of magnetic field decay 



dEAit) a^(t) dICmiB) dB{t) 



dt 



dB 



dt 



(14) 



In the reminder we focus on the conversion of energy 
of Alfven vibrations into power of its magneto-dipole 
radiation which is described by equation 



dEA[t) 
dt 



3c3 



(15) 



Let us consider torsional magneto-mechanical oscilla- 
tions which are accompanied by fluctuations of total 
magnetic moment 5xr\.[t) preserving its initial (in seis- 
mically quiescent state) direction m = m n = constant, 
i.e. a case when <5m(t)||m. The total magnetic dipole 
moment should execute oscillations with frequency uj{t) 
equal to that for magneto-mechanical vibrations of stel- 
lar matter which are described by equation for a{t). 
This means that 5m{t) and a{t) must be subjected to 
equations with similar form, namely 



5m.{t)+uj'^(t)6m.{t) = 0, 



(16) 



a{t)+u'{t)a{t)^0, u\t)^B\t)^. (17) 



It is easy to see that equations p6| and (fT7| can be 
reconciled if 



5m{t) = m a{t). 



(18) 



With account of all this, the equation is reduced 
to the equation of the magnetic field decay 



dB{t) 



^iB^it) 



dt 

T^[2^B\Q)]-\ 



1 



B(t) = 

2 2i?772 



B(0) 



V^Tt/i 



9MMc^ 



(19) 



(20) 



This shows that duration of decay, r strongly depends 
upon intensity of initial field B{0): the larger 5(0), the 
shorter decay time r. 

3.1 Magnetic-field-decay-induced lengthening of pulse 
period 

The most striking effect of magnetic field decay on 
vibration-powered radiation is the lengthening of peri- 
ods of pulsating magneto- dipole radiation, as is demon- 
strated in Fig.l. The period and its derivative are given 

by 

P(0) 



P{t) = P(0) V^ + {t/T), P{t) 



2P{t)p{ty 



P\Q) 



' 2T^l + {t/Ty 
47r2 3A/ 



S2(0) Rif 



(21) 



Equating two independent estimates of r, given by 
equations (|20| and pip , one finds that relation between 
magnitude of total magnetic moment of the star m and 
P{t) and P{t) is given by 



Ad P{t)P{t), 



~ constant. 



(22) 



For a neutron star undergoing quadrupole node-free 
torsion Alfven vibrations about axis of dipole magnetic 
moment, the mass parameter Ai (represented in terms 
of moment of inertia / = {2/h)MR?), the magneto- 
mechanical stiffness /C, frequency u are given by 

M^\l^fMRMC^lB^R^.^{^i2,) 



uj{t) _ B{t) 



27r 



2n 



7R 
3M' 



(24) 



In these equations, the star mass M and radius R serve 
as input parameters. For a neutron star with M = 
1.4 Mq and radius R = 10 km, we obtain 

iy{t) = 4.6 X IQ-^^B {M/Mq)-'^^^rI^^ , Hz (25) 
m = 3.8 X 10^^ Jp{t)P{t), Gcm^. (26) 
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These are unique to a neutron star emitting pulsed 
magneto-dipole radiation at the expense of energy of 
torsion Alfven vibrations. Taking from observations 
P{t) — I'^^it) and P{t) one can extract intensity of 
magnetic field and the magnitude of total dipole mag- 
netic moment of the star. Knowing B and m, from 
equation [201 one can obtain the time of magnetic field 
decay r. The practical usefulness of these estimates is 
that they can be used as a guide in search for finger- 
prints of vibration powered neutron stars. 

The basic conclusion of the model regarding the 
lengthening of vibration period (equal to period of elec- 
tromagnetic pulse) can be demonstrated by solution of 
equation for vibration amplitude 



forr 



HO) 



l + t/r' 



(27) 



Making us of new variable, s = 1 + t/r, the above 
equation takes the form 



sa"{s)+ (3^a{s) ^0, 13^^uj^{0)t^ 



(28) 



whose analytic solution reads (e.g., Polyanin and Zait- 
sev, 2004) 

ais) = Ji(2/3si/2) + C2Yii2l3s'/^)} (29) 

where Ji(2/3si/2) ^nd Yi(2/3.si/2) arc Besscl functions 
(Abramowitz and Stegun 1972). The arbitrary con- 
stants Ci and C2 can be eliminated from the following 
boundary conditions a(t = 0) = api ct{t = t) = 
where ag the above defined amplitude at the initial 
state of vibrations (before magnetic field decay). As 
a result, the solution of ([28| can be represented in the 




Fig. 2 Time evolution of tlie vibration amplitude, a{t) 
computed with indicated with indicated parameters C, j3 
and 77. 



a(t) = C[l + (</r)]V2 

x{M2/3[l + it/r)]'/') 
Mz{t)) 



(30) 
(t/r)]^/^)}, 



Y,{z{T)y 



c 



vYi{2p [1 
ao[Jiiz{0)) - 7jYi{z{0))]-^ 



The lengthening of vibration period is illustrated in Fig. 
2. in which the vibration amplitude is plotted as a 
function oi x = t/r with fixed value of C and different 
values of parameters 77 and /?, respectively. This figure 
shows that rj is the parameter regulating magnitude of 
amplitude: the larger rj the higher the amplitude. How- 
ever, this parameter does not affect the rate of period 
lengthening. Both, the elongation rate of vibration pe- 
riod and magnitude of vibration amplitude are highly 
sensitive to parameter /3. 

It worth emphasizing that the magnetic-field-decay- 
induced loss of vibration energy is substantially dif- 
ferent from the vibration energy dissipation caused by 
shear viscosity of matter resulting in heating of stellar 
material. The characteristic feature of this latter mech- 
anism of vibration energy conversion into the heat (i.e., 
into the energy of non-coherent electromagnetic emis- 
sion responsible for the formation of photosphere of the 
star) is that the frequency and period of vibrations are 
the same as that in the case of viscous-free vibrations 
(Bastrukov et al. 2010). However, it is no longer so 
in the case under consideration. It follows from above 
that the magnetic field decay resulting in the loss of 
total energy of Alfven vibrations of the star causes its 
vibration period to lengthen at a rate proportional to 
the rate of magnetic field decay. 

3.2 Comparison with model of rotation-energy 
powered neutron star 

All the above shows that in the model of vibration- 
energy powered magneto-dipole emission under consid- 
eration, the equation of magnetic field evolution is ob- 
tained in similar fashion as equation for the angular 
velocity fi(i) does in the standard model of rotation- 
energy powered magneto-dipole radiation 



dE 2 .2 



E 



I = -M R-" 
5 



(31) 
(32) 



One of the basic postulates of the model of rotation- 
energy powered emission is that the time evolution of 
(5m is governed by the equation 



6ni{t) = ['^{t) X ['^(t) X m]]) m = constant, (33) 
Sifi^ = m^il'*, m± = m sin 9 (34) 



Vibration powered neutron star 



5 



where 9 is angle of inclination of m to Q,{t) . The utilized 
in this model interrelation between the undisturbed to- 
tal magnetic moment m and constant magnetic field 
B frozen in the star presumes that the neutron star 
matter is in the state of permanent magnetization M 
of non-fcrromagnctic type: M ~ (3/87r)B (sec, for in- 
stance, Landau, Lifshitz, Pitaevskii 1995 [§76, problem 
2]), namely 



M dV 



1 



-BR^ 



B 



2m 



(35) 



This estimate shows that in rotation powered neutron 
star, the frozen-in the star magnetic field operates like 
a passive promoter of magneto-dipole radiation, that is, 
remain constant in the process of radiation. As a re- 
sult, the equation of energy conversion from rotation to 
radiation is reduced to equation for rotation frequency 



»(0) 



K 



_ 2m\ 
~ 37^' 

1 = 2K^\Q). 



(36) 
(37) 



Thus, in the model of vibration-energy powered emis- 
sion the elongation of pulse period (and, hence, the ori- 
gin of P) is attributed to magnetic field decay, whereas 
in the model of rotation-energy powered emission the 
lengthening of pulse period is ascribed to the slow down 
of the neutron star rotation (e.g., Manchester & Tay- 
lor 1977, Lorimer & Kramer 2004, Bisnovatyi-Kogan 
2010). 



4 Summary 

Since celebrated hypothesis of Baade-Zwicky about 
birth of neutron stars in core- collapse supernova and 
Wolticr-Ginzburg proposition that they should emerge 
as sources of super strong dipolc magnetic fields, it has 
been realized that such a field should operate as a chief 
promoter of magneto-dipole electromagnetic emission 
whose energy supply can be maintained either by en- 
ergy of rotation or by energy of Alfvcn vibrations (e.g.. 



Table 1 Basic frequency of Alfven vibrations, ua, and the 
time of magnetic field decay, r, computed with parameters 
of typical pulars and magnetars. 



M{Mq) 


i?(km) 


B(G) 


!^a(Hz) 


•^(yr) 


0.8 


20 


1012 


3.2510-3 


4.5310" 


1.0 


15 


10" 


2.5210-2 


2.9810^ 


1.1 


13 


10" 


0.22 


7.4 10^ 


1.3 


11 


10" 


1.89 


2.38 



Pacini 2008). As we mentioned, the idea that neutron 
stars can radiate like Hertzian magnetic dipole at the 
expense of energy of hydromagnetic vibrations has first 
been considered by Hoyle ct al. (1964), but without 
discussion of peculiarities of vibration-energy powered 
emission. With this in mind, one of the prime goals of 
this work was to explore some characteristic features of 
such emission. What is newly disclosed here is that nec- 
essary condition for the energy conversion from Alfven 
torsion vibrations into magneto-dipole radiation is the 
decay of internal magnetic field. The basic prediction of 
developed theory is the lengthening of periods of pul- 
sating vibration-energy powered magneto-dipole radi- 
ation which is caused by decay of internal magnetic 
field. Taking into account that both models of rotation 
powered and vibration powered pulsating emission of 
neutron stars predict lengthening of pulse period it is 
of interest to single out the potential subclass of objects 
whose pulsating radiation could be explained by vibra- 
tions, rather than rotation. In so doing, in Table 1 we 
present estimates of va and r as functions of increasing 
intensity of magnetic field for neutron stars with typi- 
cal masses and radii. It is seen that both above char- 
acteristics arc quite sensitive to M and R and, thus, 
to equation of state of stellar matter. For magnetic 
fields of typical radio pulsars, B ~ 10^2 the com- 
puted frequency va is much smaller than the detected 
frequency of pulses. This shows that considered mecha- 
nism of vibration powered pulsating radiation has noth- 
ing to do with pulsating emission of radio- pulsars. In 
the meantime, for neutron stars endowed with mag- 
netic fields B ^ IQi"' G our estimates of va fall in the 
realm of observed frequencies of high-energy pulsating 
emission of soft gamma repeaters (SGRs), anomalous 
X-ray pulsars (AXPs) and sources exhibiting similar 
features. According to common belief, these are mag- 
netars - highly magnetized neutron stars whose radia- 
tive activity is related with magnetic field decay (e.g.. 
Woods and Thompson 2006). As was emphasized, the 
decay time of magnetic field (and, hence, the duration 
of vibration powered pulsating radiation) strongly de- 
pend on the intensity of initial magnetic field of the 
star: the larger magnetic field the shorter its decay 
time. Also, the magnitude of r indicate that definite 
conclusion about secular lengthening of pulse period of 
vibration powered neutron star can be made only on 
the basis of reliable statistics which demands fairly long 
monitoring time, namely, several years if not decades. 
All these features of the vibration powered radiation 
of a neutron star are consistent with data on long- 
term monitoring activity of AXP, IE 2259.1+586 (2 
decades, 1980-2000) with using Rossi X-ray Timing Ex- 
plorer (RXTE) (Gavriil et al. 2004). Similar trend in 
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decreasing frequency of long-periodic pulsating emis- 
sion has been revealed in RXTE observations of the 
source XTE J1810-197 (Ibrahim et al. 2004). However, 
the opposite tendency (increase of pulse frequency and, 
hence, shortening of period) has been disclosed in activ- 
ity of the X-ray source CXO J164710. 2-455216 (Muno 
et al. 2006) after short time of its monitoring (May 
and July, 2005) with using of Chandra facility. The 
most plausible reason of this discrepancy (with both, 
the above observations and theoretical expectations of 
considered here models) is too short time of monitoring 
or, in other words, the lack of reliable statistics. With 
all above, we conclude that vibration powered neutron 
stars, if exist, can be revealed among the magnetars - 
AXPs/SGRs - like sources whose persistent X-ray lumi- 
nosity, lO^'' < Lx < 10^^ ergs s~^, cannot be explained 
by the loss of energy of slow rotation. 
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